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Abstract 

There are at least four versions of Hofava-Lishitz gravity in the literature. 
We consider the version without the detailed balance condition with the 
projectability condition and address one aspect of the theory: avoidance of 
caustics for constant time hyp ersurf aces. We show that there is no caustic 
with plane symmetry in the absence of matter source if A 7^ 1. If A = 1 is a 
stable IR fixed point of the renormalization group flow then A is expected 
to deviate from 1 near would-be caustics, where the extrinsic curvature 
increases and high-energy corrections become important. Therefore, the 
absence of caustics with A ^ 1 implies that caustics cannot form with this 
symmetry in the absence of matter source. We argue that inclusion of 
matter source will not change the conclusion. We also argue that caustics 
with codimension higher than one will not form because of repulsive gravity 
generated by nonlinear higher curvature terms. These arguments support 
our conjecture that there is no caustic for constant time hypersurfaces. 
Finally, we discuss implications to the recently proposed scenario of "dark 
matter as integration constant" . 



1 Introduction 



Hofava recently proposed a power-counting renormalizable theory of gravitation [H[2]. 
Since one of the most important aspects of the theory is a Lifshitz-type anisotropic 
scaling, it is often called Hofava-Lifshitz gravity. Various aspects of this theory were 
investigated [3J-[59J. 

Actually, in the literature there are at least four versions of the theory: with/without 
the detailed balance condition; and with/without the projectability condition. Among 
them, only the version without the detailed balance condition with the projectabil- 
ity condition has a potential to be theoretically consistent and cosmologically viable 
while there still are many unsolved issues. Hofava's original proposal [lj was with the 
projectability condition and with/without the detailed balance condition 0. 

The detailed balance condition restricts the form of potential in the 4-dimensional 
Lorentzian action to a specific form in terms of a 3-dimensional Euclidean theory. 
From cosmological viewpoint, this condition leads to obstacles [T2j W2\ and thus must 
be abandoned. 

On the other hand, the projectability condition stems from the fundamental sym- 
metry of the theory, i.e. the foliation-preserving diffeomorphism invariance, and 
thus must be respected. The foliation-preserving diffeomorphism consists of the 3- 
dimensional spatial diffeomorphism and the space-independent time reparametriza- 
tion. Since the lapse function is essentially the gauge degree of freedom associated 
with the time reparametrization, it is natural to restrict it to be space- independent. 
This is the projectability condition, and one can easily show that the foliation- 
preserving diffeomorphism preserves this condition. The projectability condition im- 
plies that the Hamiltonian constraint is not a local equation satisfied at each spatial 
point but an equation integrated over a whole space. Abandoning the projectability 
condition and imposing a local version of the Hamiltonian constraint would result in 
phenomenological obstacles [36] and theoretical inconsistencies [37]. Those problems 
disappear once the projectability condition is respected and if only the global Hamil- 
tonian constraint is imposed (see, for example, section 5 of [H]). Note that Hofava's 
original proposal was with the projectability condition. 

For these reasons, in the present paper we restrict our attention to the version 
without the detailed balance condition with the projectability condition. We suppose 
that the dynamical critical exponent z in the ultraviolet (UV) is equal to or larger 
than 3 since power-counting (super-)renormalizability requires it. To the best of the 

hofava put much emphasis on the detailed balance condition but considered it as just a way to 
reduce the number of independent coupling constants. 
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author's knowledge, no explicit inconsistency has been found against this version q 

This version of Hof ava-Lifshitz gravity, with general z, has an interesting cosmo- 
logical consequence even in the infrared (IR). The global Hamiltonian constraint is 
less restrictive than its local version, and allows a richer set of solutions than in gen- 
eral relativity. Actually, a component which behaves like pressureless dust emerges 
as an "integration constant" of dynamical equations and momentum constraint equa- 
tions [41j. As a result, classical solutions to the infrared limit of Hof ava-Lifshitz 
gravity can mimic solutions to general relativity plus cold dark matter. We shall 
discuss more about the "dark matter" in Sec. HI 

Also in the UV, cosmology based on Hofava-Lifshitz gravity has number of inter- 
esting properties. The anisotropic scaling with the dynamical critical exponent z = 3 
leads to a new mechanism of generating scale-invariant cosmological perturbations [S]. 
The nonlinear higher spatial curvature terms lead to regular bounce solutions in the 
early universe [HE]. It was also suggested that nonlinear higher spatial curvature 
terms might make the flatness problem milder than in general relativity [5]. Since 
these mechanisms do not rely on the detailed balance condition or a local version of 
the Hamiltonian constraint, they can be applied to the version we are considering 
now. 

Now, while cosmological consequences are interesting and it is worthwhile ex- 
ploring further, one has to be aware that there are many fundamental issues to be 
addressed in the future. First of all, renormalizability beyond power-counting has 
not been proven. Second, the renormalization group (RG) flow of various coupling 
constants has not been investigated. In particular, recovery of general relativity in 
the IR relies on the assumption (or hope) that the parameter A (see the next section) 
flows to 1 in the IR. Without knowing the condition for this behavior to be realized, 
we cannot be sure about recovery of general relativity in the IR. Third, this theory 
has not yet been intended to be a part of unified theory. Clearly, further develop- 
ments or/and embedding into a "bigger" theory (or other way around) is needed. In 
particular, since the "limit of speed" is subject to the RG flow [57J, we need a new 
idea to ensure that different species including those in the standard model of parti- 
cle physics are somehow related to each other so that their "limits of speed" agree 
with the "velocity of light" within experimental limits Fourth, an analogue of the 

2 Ref. [55 argued that the version with the projectability condition is also inconsistent, based 
on the following two claims: (i) constant-time hypersurfaces form caustics; (ii) the scalar graviton 
gets strongly coupled at all scales in Minkowski background even away from A = 1. Actually, these 
claims are not correct. See the last paragraph of subsection 13.11 and the fourth-to-the-last and 
third-to-the-last paragraphs of section [4] of the present paper. 

3 See e.g. refs. [60j [6TJ [62] for tight experimental limits on Lorentz violation. 
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Vainshtein effect [63] must be investigated for the scalar graviton in the limit A — > 1. 
Unlike general relativity, Hofava-Lifshitz gravity has not only a tensor graviton with 
two polarizations but also a scalar graviton. The dynamics of the scalar graviton and 
its fate in the IR must be investigated in details. Since the time kinetic term (together 
with gradient terms) of the scalar graviton vanishes in the A — ► 1 limit [H |2] , one has 
to take into account nonlinear interactions to see if it really decouples. 

Aside from those issues, there is another important question. In general relativity, 
flow of pressureless dust generically forms caustics. Thus, one might expect that the 
flow of "dark matter as integration constant" would also form caustics. If this were 
the case then constant time hypersurfaces would develop singularities since the flow of 
"dark matter" is orthogonal to constant time hypersurfaces. This would be a disaster 
for Hofava-Lifshitz gravity with the projectability condition since the constant time 
hypersurfaces have physical meaning in this theory. Indeed, if caustics formed then 
the extrinsic curvature would diverge. 

In this paper, we shall argue that this naive expectation is not correct and con- 
jecture that there is no caustic for constant time hypersurfaces. While the proof of 
this conjecture is beyond the scope of this paper, we shall provide supporting argu- 
ments. In particular, we shall show that there is no caustic with plane symmetry 
in the absence of matter source if A ^ 1. Since near would-be caustics the system 
enters the UV regime and A is expected to deviate from 1 via RG flow 0, this im- 
plies that caustics with codimension one do not form in the absence of matter source. 
We shall argue that inclusion of matter source will not change the conclusion. Since 
caustics with lower codimensions are more difficult to bounce than those with higher 
codimensions [65], this result provides a strong support for the conjecture. We shall 
also argue that caustics with codimension higher than one will not form because of 
repulsive gravity generated by nonlinear higher curvature terms. Finally, we shall 
discuss implications to the recently proposed scenario of "dark matter as integration 
constant" . 

The rest of this paper is organized as follows. In Sec. [2] we summarize basic 
equations in Hofava-Lifshitz gravity with the projectability condition. In Sec. [3] we 
shall argue that caustics do not form. Sec. 0] is devoted to discussion of implications 
to the "dark matter as integration constant" scenario. 

4 Here, it is assumed that the theory is renormalizable. It is also assumed that A = 1 is a stable 
IR fixed point of RG flow. In this case, by reversing the direction of the RG flow and going towards 
the UV, A should deviate from 1. 
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2 Basic equations 



The basic quantities in the theory are the spatial metric </y , the shift vector N l and the 
lapse function N. While and N l can depend on both space and time coordinates, 
iV can depend only on t. The fundamental symmetry of the theory is the invariance 
under the foliation-preserving diffeomorphism: 

f->t(i), x'^x'fax). (2.1) 

Under the infinitesimal transformation 

5t = f(t), 5x* = C(t,x) } (2.2) 

gij, N l and N transform as 

fyij = f d t9ij + 

5N* = d t (N*f) + d t C + C C N\ 

6(Ni) = d t (NJ) + 9ij d t C j + JC^N,, 

5N = d t (Nf). (2.3) 

Thus, N remains independent of spatial coordinates after transformation. In the IR, 
where dt and dx l have the same scaling dimension, it makes sense to assemble py, N l 
and N into a 4-dimensional metric in the ADM form: 

ds 2 = -N 2 dt 2 + g ij (dx i + N i dt)(dx j + N j dt). (2.4) 

The action is 

M 2 r 

I = I g + I m , I g = ^ L J dtdx 3 N^(K^K l3 -\K 2 + R + L z>1 ), (2.5) 

where 



Kij = ^( d t9ij - D i N j ~ DjNi), K = g^K ij: (2.6) 

Di is the covariant derivative compatible with g^, R is the Ricci scalar of g^, L z>x 
represents higher spatial curvature terms and I m is the matter action. Here, we have 
rescaled the time coordinate so that the coefficients of K' lj Kij and R agree. The 
cosmological constant term will be included in matter action if necessary. Note that 
not only the gravitational action I g but also the matter action I m should be invariant 
under the foliation-preserving diffeomorphism. 

By variation of the action with respect to N(t), we obtain the Hamiltonian con- 
straint 

H g± + H m± = 0, (2.7) 
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H g± = -j± = J dx^H g± , 7i a± = -^(K i *p ij -R-L, >1 ), 

H m± = - 5 j^ = J dx 3 Jg Tl Ti = T^n\ (2.8) 
Here, Pij and n M are defined as 

Pij = K i:j - \Kg i:j , (2.9) 

and 

n^ = -Ndt, n»d li = -(d t -N i d i ). (2.10) 
Variation with respect to N l (t,x) leads to the momentum constraint 

Hgi + H m i = 0, (2-H) 

n " s ~iiw = - M " D '^ 

As in general relativity, the gravitational action can be written as the sum of 
kinetic terms and constraints up to boundary terms: 

I g = J dtdx 3 [n ij d t gij - N i H gt ] - J dtNH g± + (boundary terms), (2.13) 

where 7r y is momentum conjugate to gij given by 



7T 



A T 

V = —i - Ml., fnrij rjj = n ik J l 



m(j , Mpiy/9P t3 , V > ~ .'/'V/'/./- (2.14) 

The Hamiltonian corresponding to the time t is the sum of constraints and boundary 
terms as 

H g [dt] = NH g± + J daPWHgi + (boundary terms). (2.15) 
Finally, by variation with respect to gij(t,x), we obtain dynamical equation 

£gij ~l" £mij = 0, (2.16) 

2 SL 



£gij — 9ik9j 



l 9 
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The explicit expression for 8„ij is given by 



e, 



hi.) 



M 2 Pl 



--(d t - N k D k ) PlJ - K Pij + 2K k p kj 



+-^(PikDjN k +p jk DiN k ) + -9ijK kl p kl - Gij 



z>liji 



(2.18) 



where £ z> uj is the contribution from L z> \ and is Einstein tensor of g^. 

The invariance of I a under the infinitesimal transformation (12.31) leads to the 
following conservation equations, where a represents g or m. 



Nd t H a± + / dx 



N i d t (^gH ai ) + -N^gE%d t g ij \ = 0, 
^(dt-NW^Hai + KHai-^HajDiW -D j £ aij = 0. 



(2.19) 



3 Absence of caustics 



The vector n 11 defined in (12.101) has unit norm and is orthogonal to constant time 
hypersurfaces. Since the lapse function iV depends only on time, follows the 
geodesic equation. 

n^V M n, = n»V u n^ = -9(n%) = 0. (3.1) 

In general relativity, a congruence of geodesies would generically form caustics. As 
an example, let us consider a congruence of geodesies orthogonal to the hypersurface 
t = T(x) in Minkowski spacetime ds 2 = —dt 2 + dx 2 + dy 2 + dz 2 . By introducing a new 
time coordinate r as the proper time along each geodesic and a new spatial coordinate 
X as the value of x at the intersection of each geodesic with the hypersurface t = T(x), 
the Minkowski metric is rewritten as 



ds 2 = -dr 2 + a 2 (r, X)dX 2 + dy 2 + dz 2 , a(r, X) = yjl - [T'(X)f 



T 



Tc(X) 



where 



Tc(X) 



[1 - T'(X) 2 f 2 



(3.2) 



(3.3) 



-T"(X) 

The metric component a 2 (r,X) vanishes at finite proper time r = t c (X) and, thus, 
the congruence of geodesies form caustics. In general relativity with a certain energy 
condition, it is easy to show that a congruence of geodesies forms caustics in more 
general situations essentially because gravity is attractive. 

On the other hand, for the vector n M in the Hofava-Lifshitz gravity, higher curva- 
ture terms become important near (would-be) caustics and provide negative effective 
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energy and repulsive gravity. Also, A should deviate from 1 by RG flow (see footnote 
Hj) and, thus, the kinetic terms also contribute differently from general relativity. 

3. 1 Codimension higher than one 

For codimension higher than one, the spatial curvature of the constant time hypersur- 
face increases near (would-be) caustics. The system enters the non-relativistic regime 
and higher spatial curvature terms become important. Among them, highest order 
terms (e.g. curvature cubic terms in the case of z — 3) generate the strongest restoring 
force. As in some early universe models [UEJE], we expect that the (would-be) caus- 
tics should bounce at short distance scales if codimension is higher than one. Note 
that this is not because of deviation from geodesies | bnt becanse of repnlsive gravity 
at short distances. In general relativity, congruence of geodesies would form caustics 
because gravity is attractive. On the other hand, for the vector in Hofava-Lifshitz 
gravity, higher curvature terms become important near the (would-be) caustics and 
provide repulsive gravity and thus bounce. 

For an odd z, the sign of the highest nonlinear spatial curvature term can in 
principle change. Nonetheless, as far as a contracting region has a finite volume, 
the leading (would-be) divergence in the spatial curvature at late time is expected 
to be positive. For this reason, even with an odd z, we expect would-be caustics 
to bounce eventually. On the other hand, numerical confirmation of this kind of 
behavior probably requires a rather wide dynamic range since we have to wait until the 
finiteness of the contracting region becomes important. Thus, for numerical purposes 
it is probably easier to consider a large enough, even z. We hope to perform numerical 
analysis of the bouncing behavior in the future. 

One might worry about the fact that the anisotropic scaling in the UV leads 
to the scaling oc l/a z+3 for radiation energy density [20]. For z — 3, this might 
cause difficulties for bouncing cosmology with FRW (Friedmann-Robertson- Walker) 
spacetime since the radiation energy density scales in the same way as the z = 3 
higher curvature terms and has the opposite sign. On the other hand, in the present 
situation, the contracting region has just a finite volume and thus radiation does not 
have to be comoving with the vector n M . Indeed, while the "dark matter" is pressure- 
less and the vector n M follows geodesies, the large pressure (P ra d = ( z /3)p ra d in the 
UV) acts as extra repulsive force for radiation and prevents radiation from following 

5 In the case of ghost condensate [64], the derivative of the scalar field responsible for the con- 
densate deviates from geodesies because of higher derivative terms [65]. On the other hand, in 
Hofava-Lifshitz gravity the vector always satisfies the geodesic equation (|3.ip . 
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geodesies during inhomogeneous contraction. Thus, radiation can easily diffuse and 
grows much more slowly than the z = 3 curvature terms. For this reason, for z = 3, 
radiation does not cause difficulties. For z > 3, highest nonlinear curvature terms 
(oc l/a 2z ) grows faster than radiation energy density even in FRW spacetime and, 
thus, radiation does not prevent the highest curvature terms from acting as restoring 
forces. 

Therefore, unless the whole universe contracts, (would-be) caustics with codimen- 
sion higher than one should bounce due to higher curvature terms. If the universe is 
completely homogeneous and flat, i.e. if the universe is a flat FRW spacetime, then 
higher spatial curvature terms vanish classically. Hence, a contracting flat FRW does 
not bounce classically. However, quantum mechanically, there must be fluctuations 
and nonlinear higher spatial curvature terms must acquire non-vanishing expectation 
values. Those fluctuations grow as the universe contracts. Therefore, if z is large 
enough and perhaps if z is even (so that the highest nonlinear spatial curvature terms 
always act as a restoring force) then a contracting flat FRW universe might also 
bounce after all. Further investigation of this issue is worthwhile. 

Note that the bounce is provided by nonlinear terms. Therefore, if we analyzed 
behaviors of the vector n M and constant time hypersurfaces without including those 
nonlinear terms then we would not be able to see the bounce and would instead see 
caustics forming. It is likely that this is closely related to instabilities of linear pertur- 
bation found in |38j. It is also important to include backreactions of the higher spatial 
curvature terms to the geometry since, as already stated, the bounce is not due to 
deviation from geodesies but due to repulsive gravity at short distances. If we did not 
take into account backreactions to the geometry, one would simply conclude forma- 
tion of caustics [55J. Without taking into account nonlinear terms and backreactions 
to the geometry, we would never be able to describe the system properly. 

3.2 Codimension one 

On the other hand, for codimension one, higher spatial curvature terms do not help. 
This is in accord with the observation in [65J that caustics with lower codimensions 
are more difficult to bounce than those with higher codimensions. In order to see it, 
let us consider the following ansatz with plane symmetry. 

N = a(t), N l di = (3(t, x)d x , g ij dx i dx j = a(t, xfdx 2 + dy 2 + dz 2 . (3.4) 

We have the foliation-preserving diffeomorphism 

t—*t(t), x— > x(t, x), y — > y, z — > z. (3.5) 
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By using this symmetry, we can set a = 1 and a — 1. The ansatz is now reduced to 



N = 1, APd; = /3{t, x)d x , g ij dx i dx j = dx 2 + dy 2 + dz 2 . (3.6) 

We still have the residual symmetry 

t— >t + t , x^x + x (t), (3.7) 

where to is a constant and xq is a function of t. It is now evident that higher spatial 
curvature terms do not help since the spatial metric is flat. 

Actually, what prevents caustics with codimension one from forming is deviation 
of A from 1. In order to recover general relativity in the IR, A = 1 must be an IR 
fixed point of RG flow. However, near (would-be) caustics, the system enters the UV 
regime and A should deviate from 1 by RG flow (see footnote H]). In the following, we 
shall show that there is no caustics with codimension one if A ^ 1. 

Ideally speaking, it is appropriate to take into account fully quantum mechanical 
effects to analyze behavior of the system near would-be caustics. However, at this 
stage where renormalizability (beyond power-counting) and RG flow of the theory 
are not yet understood, it is not easy to perform fully quantum mechanical analysis. 
For this reason, in this paper, we partially take into account quantum effects by 
including the scale-dependence of coupling constants (especially, deviation of A from 
1) in classical equations of motion. 

As already argued in the previous subsection, radiation energy density does not 
become significant near (would-be) caustics simply because radiation does not have 
to be comoving with the vector and thus can diffuse. The same applies to other 
forms of matter fields. For this reason, in the following we ignore matter fields and 
set T, 



0. 



For the ansatz (13.61) . it is easy to show that 



n 



s, 



gxx 



£, 



-(A - l)M 2 Pl /3" 
(A - 1)M 



mm 



£, 



gzz 



XM 2 Pl 



-f3' + f 3(3" + -(/3r 
-P> + {3(3"+ l -{(3') 2 



M 2 



(3.8) 



where an overdot and a prime denote derivatives with respect to t and x, respectively. 
With A 7^ 1, the general solution to the set of equations Ti gx = £ gxx = £ gyy = £ gzz = 
is 

= #,(*), (3.9) 
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where (3o(t) is an arbitrary function of t 0. Thus, by using the residual symmetry 
(13. Tp we can set 

P = 0. (3.10) 
This shows that there is no caustics with codimension one. 

4 Dark matter as integration constant 

As already mentioned, in Hofava-Lifshitz gravity there is no local Hamiltonian con- 
straint. This might cause worries since in general relativity the local Hamiltonian 
constraint is nothing but the Poisson equation, which is one of the most important 
equations for gravity. Actually, in the following we shall see that the absence of lo- 
cal Hamiltonian constraint leads to an interesting consequence. A component which 
behaves like dark matter emerges after solving the system of equations when we in- 
terpret general solutions. We will see that the Poisson equation with "dark matter" 
built-in is satisfied by the solutions. 

Now, following ref. [41J, let us define energy density p d of "dark matter as inte- 
gration constant" by 

M 2 

p d = —^{l\' J i>,j -R- L Z>1 ) - T^n". (4.1) 

Consistency with (12.191) requires that 

d t J dx 3 ^ Pd = 0, (4.2) 

but this is automatically satisfied because of the Hamiltonian constraint. Equations 
of motion in Hofava-Lifshitz gravity are summarized as 

M 2 Pl G^ v = T^ + Pd nV, (4.3) 
dx 3 ^p d = 0, (4.4) 



where 



dY Jd 



+ n 



V 



-£ ■ o ik a jl 




6 With A = 1, wc would obtain (3 = flm = /3o(t) + x/(t c — t), where Po(t) is an arbitrary function 
of t and t c is a constant. By using the residual symmetry (|3 . 7[) we can set /3jr = x/(t c — t). This 
would diverge as t — > t c and thus represents a would-be caustics. However, as already stated, near 
a would-be caustics, the system enters the UV regime and A should deviate from 1 by RG flow (see 
footnoted]). Thus, this solution is invalidated before actually reaching t = t c . 
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By taking divergence of (14. 31) . we obtain 



(d± Pd + K Pd )n, = -V U (T^ - M 2 Pl G$), (4.6) 

where d± = n^d^. Because of the spatial diffeomorphism invariance, the right hand 
side is proportional to n^. Thus, this equation has only one non- vanishing compo- 
nent. By contracting with n^, we obtain the (non-)conservation equation of "dark 
matter" [31] : 

d ±Pd + Kp d = n*V v [T^-M* l G ( $]. (4.7) 

In the IR limit with A — > 1, G^^ u reduces to the 4-dimensional Einstein tensor G^^ u 
and eq. (14. 3p reduces to the Einstein equation with "dark matter" 

M 2 piG {4)„v = T „u + pdn » n v. (4. 8 ) 

From this equation, one can obtain the Poisson equation with "dark matter" built-in. 
In the same limit, (14. 7p reduces to the conservation equation, 

d ±Pd + K Pd = 0. (4.9) 

Note that p d can be positive everywhere in our patch of the universe. In a homo- 
geneous spacetime such as the FRW spacetime, the global Hamiltonian constraint is 
as good as local one since all spatial points are equivalent. However, in inhomoge- 
neous spacetimes there are drastic differences. If the whole universe is much larger 
than the present Hubble volume then it is natural to expect that the universe far 
beyond the present Hubble horizon is different from our patch of the universe inside 
the horizon. In this case, the global Hamiltonian constraint (14.41) does not restrict 
the universe inside the horizon. Even if we approximate our patch of the universe 
inside the present horizon by a FRW spacetime, the whole universe can include in- 
homogeneities of super-horizon scales and, thus, the global Hamiltonian constraint 
does not restrict the FRW spacetime which just approximates the behavior inside the 
horizon. Therefore, p d can be positive everywhere in our patch of the universe inside 
the present Hubble horizon. 

In the UV epoch, fluctuations of matter fields as well as metric fluctuations act 
as the source term in (14. 7p . Since such fluctuations include modes with various wave- 
lengths, fluctuations of "dark matter" are generated with various wavelengths, in- 
cluding those far longer than the size corresponding to the present Hubble horizon. 
Therefore, there must certainly be large enough regions with positive p d . In principle 
it should be possible to predict a typical amplitude of p d , once a model of the early 
universe is specified in the context of Hofava-Lifshitz gravity. 
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Note also that we do not have to promote the "dark matter as integration con- 
stant" to an independent dynamical field as far as the scalar graviton, the tensor 
graviton and matter fields are considered as independent dynamical fields in the initial 
value formulation. The initial value formulation of Hofava-Lifshitz gravity consists of 
dynamical equation (12.161) . global Hamiltonian constraint (12.71) . local momentum con- 
straint (12.111) and gauge conditions. Of course, the constraints are preserved by the 
dynamical equation. In this language, the "dark matter" emerges only after solving 
the system of equations when we try to interpret a solution. 

In order to describe the scalar graviton, the commonly used method called Stiickelberg 
formalism does not seem to be useful here. If we adopted it to construct an effective 
field theory of the scalar graviton then, unlike ghost condensate [64] . the foliation- 
preserving diffeomorphism would forbid h^ and thus n 2 . For this reason, this de- 
scription would not include a healthy kinetic term even with A ^ 1 and would get 
strongly coupled at all scales in Minkowski background [55] . In non- vanishing "dark 
matter" backgrounds, the strong-coupling scale of the Stiickelberg field becomes fi- 
nite but is still as low as p 1 / 4 even with A ^ 1 [55J. This indicates breakdown of this 
description but does not imply inconsistency of the underlining UV theory. The phys- 
ical reason for this is that the Stiickelberg field has vanishing overlap with the scalar 
graviton. Indeed, if we adopt for example the gauge used in p] (without introducing 
a Stiickelberg field) then the scalar graviton has a finite kinetic term away from A = 1 
and does not exhibit the strong coupling mentioned above for the Stiickelberg field. 
(If we introduce the Stiickelberg field then we should use a gauge-invariant variable 
representing the scalar graviton and solve constraint equations. After all, we should 
be able to obtain a finite kinetic term away from A = 1 in Minkowski background.) Of 
course, even in this description, we have to take into account nonlinear interactions 
carefully when we take the limit A — > 1 to see if there is an analogue of the Vainshtein 
effect [63]. We hope to come back to the issue of Vainshtein-like effect in the near 
future. 

If we consider a group of many microscopic lumps of "dark matter" then collisions 
and bounces among those lumps may accumulate to generate non-trivial effects in 
macroscopic scales. This is exactly the spirit of renormalization group. As far as 
gravity at astrophysical scales is concerned, this kind of collective behavior of small 
lumps of "dark matter" might mimic behavior of a cluster of particles with velocity 
dispersion. If this is the case then microscopic lumps of "dark matter" play the role of 
particles in usual dark matter models. It is certainly interesting to see what happens 
if a group of astrophysically large number of such microscopic lumps of "dark matter" 
collides with another group. Clearly, detailed investigation is necessary to understand 
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rich dynamics of "dark matter" from microscopic to macroscopic scales. 

The foliation-preserving diffeomorphism invariance is the fundamental symmetry 
of the theory. Thus, the whole system including all matter fields must respect this 
symmetry. In particular, the matter action must be invariant under the 3-dimensional 
spatial diffeomorphism. This fact, combined with the orthogonality of the flow of 
"dark matter" to the constant time hypersurface, means that the dispersion relation 
for each matter field, 

- 2 = ^^ + --- + ^ 2 + m 2 , (4.10) 

is defined in the rest frame of "dark matter" . Thus, in the region where matter fields 
move relative to "dark matter" with large relative velocities, higher order terms in the 
dispersion relation of matter fields can become important. It is certainly interesting 
to investigate astrophysical implications of such effects. 

Note added 

The issue of caustics in Hofava-Lifshitz gravity with the projectability condition was 
discussed in [55 J. Following recommendation by an anonymous referee, here we would 
like to comment on it. Ref. [55] indeed has three statements about Hofava-Lifshitz 
gravity with the projectability condition: (i) "dark matter" forms caustics; (ii) "dark 
matter" is described by ghost condensate |64j ; (iii) the scalar sector gets strongly 
coupled at the scale A ~ p 1 / 4 even with A ^ 1. Actually, these three comments are 
not correct for the following reasons, (i) They did not take into account repulsive 
gravity due to nonlinear higher curvature terms as mentioned in Sec. E]of the present 
paper, (ii) Ghost condensate and Hofava-Lifshitz gravity have different symmetries 
as mentioned in the third-to-the-last paragraph of Sec. 0] of the present paper, (iii) 
The strong coupling away from A = 1 found in [55] indicates sickness of their descrip- 
tion, i.e. the way dynamical degrees of freedom are identified, but does not imply 
inconsistency of the underlining UV theory. For this point, see the fourth-to-the-last 
and the third-to-the-last paragraphs of Sec. 0] of the present paper. 
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